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Abstract 

We present three-loop formulas for the additive renormalization constant 
A(u, e) and associated renormalization group function B(u) for the specific 
heat of the O(n) symmetric tp 4 model. Using this result, we obtain also the 
amplitude function above Tq within the minimally renormalized theory at 
fixed d = 3. At the fixed point, the three-loop correction to B(u) turns out 
to be small (about 3% for n = 2). We note that a correction of this size may 
become important at the level of accuracy expected in future experiments. 

PACS: 05.70.Jk, 75.40.Cx, 67.40.Kh, ll.10.Gh 

Field-theoretic renormalization group (RG) calculations based on (Borel) resummations 
of several orders of perturbation theory have yielded accurate predictions for the critical 
exponents [HJ and for many of the universal amplitude ratios || of 0(n) symmetric systems. 
For n > 1, however, similar predictions for amplitude ratios involving quantities denned 
below Tc are not available since the relevant perturbation series have not yet been extended 
to sufficiently high order for resummations to be effective. It has been pointed out recently 
0, that such higher order calculations for the amplitude functions of the specific heat and 
supernuid density would be needed for a fully quantitative test of unversality along the 
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A-line of 4 He. Another quantity which enters the formulas for the amplitude ratios and 
which has, so far, been computed only in low order is the renormalization group function 
B{u) associated with the additive part of the renormalization of the specific heat. Like the 
amplitude functions, this function has additional relevance for the analysis of experimental 
data in the nonasymptotic region 0-0]. 

On the basis of specific heat measurements taken in earth orbit, Lipa et al f| have 
shown that the rounding of data near the A-transition of 4 He due to gravity-induced pressure 
gradients |§, can be avoided for reduced temperatures as small as t ~ 10~ 9 . Previously, these 
effects restricted the range of useful data to temperatures 10~ 6 t ^ 10~ 2 implying the need 
for theoretical constraints in the analysis ||10|| . An unconstrained fit to the data in Ref. || 
yielded the critical exponent value a = —0.01285(38) and an estimate A + /A" = 1.054(1) for 
the ratio of leading amplitudes. The uncertainty in this value of a is smaller than that of the 
"best" RG prediction [[LlJ] a = —0.016(6) by about an order of magnitude — a clear call for 



greater theoretical accuracy. Further experiments in reduced gravity have been planned fll2 |. 

The purpose of this note is to examine the three- loop approximation to B(u) from the 
standpoint of the minimally renormalized 4 theory in fixed dimension P,|T3"|-|Ttj|] . This func- 
tion is expected to deviate from its leading order approximation B(u) ~ n/2 by a small 
amount of 0(r]) [^ |14| , |T7| , where r] is the exponent describing the decay of spatial correla- 
tions at T = Tq- While the role played by B(u) is expected to be a minor one relative to the 
exponents and amplitude functions , it is nevertheless of interest to know whether 

a correction of this size would become significant at the higher level of accuracy expected in 
reduced gravity experiments. 

The renormalization of the specific heat within the minimal subtraction scheme at fixed 



dimension has been described in detail in Ref. |15|]. For definitions and notation, see also 
Refs. P,|l3|-|T6|1 . In three-loop order, the additive renormalization A(u,e) and RG function 
B{u) are given by 
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A(u, e) = -In- - 8n(n + 2)- + a— + 0(« 3 ), 



r e 2 e 3 



a = -32n(n + 2) 
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which we have obtained within the "massless" theory (that is, for k ^ and T = Tq) since 
the pole terms are more readily evaluated there. The relevant vacuum diagrams with two 
(f) 2 insertions are shown in Fig. [IJ their contributions near d = 4 are 
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To obtain J^, we have used Eq. (2.20) of Ref. ||18|| . The geometric factor A d = 2 "ir 
d/2)/(d — 2) is left unexpanded ||. In three dimensions, our formula for the amplitude 
function above Tc reads 



F + (l,u, 3) 
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For n = 2, the 0(w 2 ) term in Eq. (|2|) is roughly 3% of the leading term and, since r\ ~ 0.04 
|T3| |, is consistent with the 0(7/) estimate of Ref. |T7] for the net contribution of all higher 
order terms. It has been suggested j5|, [T^JT7|| that the terms beyond leading order should 
contribute less than 1% to the function B{u) and yet, although this contribution is expected 
to be small, it is not at all clear that it should be so small. One should bear in mind that 
low order perturbative expressions, such as Eq. (E|), cannot by themselves be regarded as 



reliable in a purely quantitative sense and that it is usually difficult to anticipate which (low) 
order of perturbation theory will provide the "best" approximation in any given situation. 
Indeed, this is the motivation behind the resummations of higher order series that have so 
far yielded accurate predictions for the exponents and amplitude ratios. 

With the above caveat, therefore, let us consider the 0(w 2 ) term in Eq. (g) to be 0(r/) 
and examine its effect on the amplitude ratios JT5 ] 
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where Q±, F_ and G are the amplitude functions for the correlation lengths, the specific 
heat below Tc and the superfluid density, respectively; the asterisk denotes fixed point 
values. We make use of the Borel summation results given in Refs. [|l^,[14],|nj and of the 
relations Q* + = 2vPt and Q*_ = 3 — 2Q1 where P + is the amplitude function for the quantity 



(dr /d£- 
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13]. In the absence of Borel results for n > 1 below Tc, we use the most reliable 



low order approximations for F_ and G, which turn out to be given already in one-loop order 
U. We also set u* = 0.0405, a = 0.11 (n = 1) and u* = 0.0362, a = -0.013 (n = 2) and fix 
v according to v — (2 — a)/3. 

The values of A + / A~ and Rj given in Table | illustrate the size and direction of the 
effect of including the 0(u 2 ) term in Eq. (||]) when all other quantities in Eqs. @ and C [3~0|) 
are kept fixed. Also shown in the table is the value for A + j A~ (n — 1) obtained by Bagnuls 
et al ||19|| , who use a different renormalization scheme, and the experimental values for 4 He 
for A + /A~ and Rj obtained by Lipa et al and by Singsaas and Ahlers respectively. 
In each case, the effect of the 0(w 2 ) term in Eq. (0) is comparable to the uncertainties given 
by the authors of Refs. [§,|1^,|1^]. Since the possibility of the exact value of B* differing from 
the leading term by ~ 3% cannot be ruled out and since the experimental uncertainties are 
expected to be substantially smaller in the future, it seems that a higher order calculation, 
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of the kind indicated by the analysis of Ref. for amplitude functions below Tc, may be 
needed for B(u). This conclusion, of course, presumes the future availability of improved 
estimates for the critical exponent a and for the amplitude functions below Tc- 

It may be argued that an additive renormalization is unnecesary if the specific heat is 
represented in terms of its temperature derivative dC^/dt. In this connection, we recall the 
relations fll5| 

f) F 

8AfP ± fS fi) = (e - 2( r )F ± + 4B- f3 u ^ } (11) 

where f±'°^ are the amplitude functions for dC ± /dt. These formulas, for example, enable 
the amplitude ratios to be expressed in terms of /± 3 '°' ) . In that case, resummation results for 
B{u) would be useful for an internal check of the theory However, in view of the unusually 
large uncertainty associated with the present high order result for /j 3,0 ^ [12% compared to 



< 



1% typically for other amplitude functions []14|,|16 |], the representation based on additive 



renormalization may well turn out to be the more reliable in quantitative applications. A 
higher order calculation of B[u) would be needed to answer this question. 

Finally, we note that in Ref. ||14|| , the higher order coefficients in the perturbation series 



for F + were approximated by use of Eq. (|TTD with B(u) ~ n/2. This procedure neglects the 
leading poles of A(u, e) beyond two-loop order [for example, the term ~ |e 2 in the square 
brackets of Eq. ([[])]. Using Eqs. (fj) and flTT|), we find that the resummation results for F* 
are shifted by about 2%. However, since F + enters the formulas for the amplitude ratios and 
for the analysis of experimental data only in the combination aF + , the effect here is entirely 
negligible. 

In summary, we have computed, within the framework of the minimally renormalized </> 4 
theory at fixed dimension d = 3, the three- loop correction to the additive renormalization of 
the specific heat, Eq. (fj), for systems with 0(n) symmetry. We have used this to determine 
the corresponding RG function B{u) in Eq. (H), and amplitude function F + above T c in 
Eq. (H|). While the neglect (within the present scheme) of the leading additive poles in the 



specific heat beyond two-loop order is justified in analyses based on low order perturbation 
theory, these poles may lead to a small systematic effect at the level of accuracy expected in 
future experiments f8|.|T2|j. 
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FIGURES 



TABLES 
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TABLE I. Values for the amplitude ratios 
obtained using B(u) ~ n/2 [1&2 loop] and 
Eq. (0) [3 loop]. All other quantities in Eqs. (|9|) 
and (|l(j) are held fixed as described in the text. 
The other theoretical and experimental values 
are included to illustrate the current level of un- 
certainty for these quantities. 



A+/A- 



n 



n 



0.527 
0.536 
0.541(14) a 

1.056 
1.054 
1.054(1) & 



[1 & 2 loop] 
[3 loop] 



[1 & 2 loop] 
[3 loop] 



a Bagnuls et al |0| 
b Lipa et al [g] 

and Ahlers fTTOl 



(A-y/ 3 /k 



0.831 
0.840 
0.85— 0.86 c 



FIG. 1. Diagrams of the massless theory 
contributing to the specific heat. Their expan- 
sions near d = 4 are given in Eqs. 
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